Nonlinear excitation of acoustic modes by large amplitude Alfven waves in a 

laboratory plasma 
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The nonlinear three-wave interaction process at the heart of the parametric decay process is 
studied by launching counter-propagating Alfven waves from antennas placed at either end of the 
Large Plasma Device (LAPD). A resonance in the beat wave response produced by the two launched 
Alfven waves is observed and is identified as a damped ion acoustic inode based on the measured 
dispersion relation. Other properties of the interaction including the spatial profile of the beat mode 
and response amplitude are also consistent with theoretical predictions for a three-wave interaction 
driven by a non-linear pondermotive force. 

PACS numbers: 52.35.Mw, 52.35.Bj 



Alfven waves, a fundamental mode of magnetized plas- 
mas, are ubiquitous in lab and space. While the linear 
behavior of these waves has been extensively studied [Ti- 
ll], non- linear effects are important in many real systems, 
including the solar wind and solar corona. In particu- 
lar, a parametric decay process in which a large ampli- 
tude Alfven wave decays into an ion acoustic wave and 
backward propagating Alfven wave may play a key role 
in establishing the spectrum of solar wind turbulence [6] . 
Ion acoustic waves have been observed in the heliosphere, 
but their origin and role have not yet been determined 
[7]. Such waves produced by parametric decay in the 
corona could contribute to coronal heating [5]. Para- 
metric decay has also been suggested as an intermediate 
instability mediating the observed turbulent cascade of 
Alfven waves to small spatial scales [SJ [S] ■ 

In this letter, the first laboratory observations of 
the Alfven-acoustic mode coupling at the heart of the 
parametric decay instability are presented. Counter- 
propagating shear Alfven waves are launched from an- 
tennas and allowed to interact nonlinearly. As the beat 
frequency between these two launched waves is varied be- 
tween discharges, a resonant response is observed when 
frequency and wave number matching is satisfied for cou- 
pling to an ion acoustic mode. Other features of the in- 
teraction including the beat mode spatial structure and 
response amplitude match predictions based on a three- 
wave interaction driven by a non-linear pondermotive 
force. 

Although these results represent a beat wave process 
rather than an instability, the reported evidence of a 
three-wave interaction may be used to validate simple 
theoretical predictions and aid in comparison with space 
measurements. To date, there has been an abundance 
of theoretical work [TDHIS]j but very little direct exper- 
imental observation of parametric decay. Observations 
by Spangler et al. [T7] in the ion foreshock region up- 
stream of the bow shock in the Earth's magnetosphere 
indicate the presence of large amplitude Alfven waves 
as well as density fluctuations with no magnetic spectral 



FIG. 1. Experimental setup in the LAPD plasma column. 
Alfven wave antennas shown at either end of the device launch 
the counter-propagating Alfven waves to be examined in this 
study. Magnetic probes B4 and B7 and Langmuir probes LI 
and L2 placed between the two antennas in the plasma column 
are used to diagnose the interaction; various z positions are 
used for the data in the paper. 



component; the latter are presumed to be acoustic modes 
resulting from parametric decay, but detailed measure- 
ments are limited. 

The Large Plasma Device (LAPD) at UCLA is an ideal 
environment for experiments diagnosing nonlinear Alfven 
wave interactions. The LAPD is a cylindrical vessel ca- 
pable of producing a 16.5 m long, quiescent, magnetized 
plasma column for wave studies. The BaO cathode dis- 
charge lasts for ~ 10 ms, including a several ms current 
flattop. Typical plasma parameters for the present study 
are n^ - lO^^ /cm^, Te - 5 eV, and Bq ~ 400 - 900 G 
(/3 ^ 1) with a fill gas of helium or hydrogen. Exten- 
sive prior work has focused on the properties of linear 
Alfven waves [5J [15^50] . Studies of the nonlinear proper- 
ties of Alfven waves have also been performed on LAPD; 
in these experiments, two co-propagating Alfven waves 
non-resonantly drive a quasi-mode in the plasma |21j or 
resonantly drive drift wave instabilities [22] . 

For the present set of experiments, loop antennas 
placed at either end of the LAPD, shown in Fig. [T] 
launch linearly polarized, counter-propagating Alfven 
waves with amplitudes of dB ~ 1 G during the discharge- 
current-flattop period of the LAPD discharge. It should 
be noted that parametric decay of a single Alfven wave is 
not observed in these experiments; consistent with this. 
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FIG. 2. Beat interaction between two counter-propagating 
Alfven waves in helium plasma showing the ring up and 
ring down associated with coupling to a natural mode of the 
plasma. Two Alfven wave antennas at frequencies of 230 and 
216 kHz are on between 8 and 10 ms. Bo = 750 G. Panel (A) 
shows the filtered ion saturation current signal between 8 kHz 
and 20 kHz from a Langmuir probe at 2 = 3.20 m. Panel (B) 
shows a singal from a nearby magnetic probe at z = 2.24 m. 
The frequency spectrum over the 2 ms window the antenna is 
on is displayed in Panel (C) for the ion saturation current sig- 
nal (upper solid trace) and the magnetic singal (lower dashed 
trace). Panel (D) shows the beat amplitude as a function of 
the product of the magnetic field amplitudes at the driving 
frequencies. All probes are at x — 0, y — 0. 



the experimental value of 6B/Bf) < 2 x 10^"^ gives a 
growth rate (TU] that is comparable to an Alfven wave 
transit time through the entire plasma column for rel- 
evant experimental parameters. Instead, antennas di- 
rectly launch both the "pump" and "daughter" Alfven 
waves at similar amplitudes. In the plasma column be- 
tween the antennas, magnetic probes detect the magnetic 
field signatures of the launched modes while Langmuir 
probes are used to detect signatures of a density response 
at the beat frequency. Each probe is mounted on an auto- 
mated positioning system that may be used to construct 
a 2-D profile in the x-y plane averaged across multiple 
discharges. 

A clear nonlinear response at the beat frequency is ob- 
served in these experiments, as shown in Fig. [2] When 
the two Alfven wave antennas are turned on between 
8 ms and 10 ms in this helium discharge, a beat wave 
at the difference frequency of 14 kHz is observed both 
in the filtered ion saturation current trace displayed in 
Panel (A) and the full frequency spectrum shown in 
Panel (C). This signal will be shown to have many prop- 
erties consistent with an ion acoustic mode produced by 
a three- wave matching process. The beat amplitude of 
75 /iA represents '^ 3.5% of the measured mean ion satu- 
ration current. After the last of the magnetic signatures 



FIG. 3. Beat amplitude as a function of beat frequency A/ 
showing a resonant response in helium plasma with back- 
ground Bo = 750 G. The Alfven wave antenna on the cath- 
ode end produces a fixed frequency 230 kHz wave while the 
frequency of the wave produced by the end mesh antenna 
is scanned from 205 to 230 kHz between discharges. The 
dashed curve shows the beat amplitude SIsat as a percent 
of Isat observed as a function of A/ for a Langmuir probe 
at z — 6.07 m. The amplitude is normalized to the zero fre- 
quency component; the shaded errorbar represents the level of 
background fluctuations. The thin dash-dot trace is the pre- 
dicted density response |pi/po| based on Eq. fusing A^ — 0.14, 
/3 = 4 X 10"" and b±/Bo = 1/750. The thin solid trace is the 
equivlent prediction for A'^ = 0.36. The flgure inset to the 
right shows the N — 0.36 trace reseated to the data to em- 
phasise that this A'' is the best fit for the width of the peak. 



from the Alfven waves pass by a fixed magnetic probe 
at t = 10.015 ms, the amplitude of the beat wave does 
not immediately drop to zero, indicating that coupling 
to a normal mode of the plasma has occurred. The ring- 
down time of the driven wave is ^ 85 /is, comparable to 
an ion-neutral collision time of ^ 100 fis for these param- 
eters. When these experiments are repeated in hydrogen 
plasmas, the ring down time is shorter; consistent with 
this, the ion- neutral collision frequency for the chosen 
parameters is larger in hydrogen than it is in helium }23j . 

The beat amplitude is expected to be largest when 
three-wave coupling most efficiently excites a normal 
mode of the plasma. The experimental strategy to test 
this prediction is as follows: the launch frequency of the 
cathode side antenna is held fixed while the launch fre- 
quency of the end mesh side antenna is varied between 
discharges. The plasma response at the beat frequency is 
then examined in each discharge to find the difference fre- 
quency that best couples to an acoustic beat mode. The 
results of this scan for helium plasmas with Bq — 750 G 
are represented by the dashed line in Fig. [3j This curve, 
representing the beat wave amplitude with both anten- 
nas on plotted as a function of the difference frequency. 



peaks at a frequency of around 13 kHz. This frequency 
at which three-wave matching relations are best satisfied 
to excite a normal mode of the plasma is defined as the 
resonance frequency. 

A calculation based on the ion acoustic and Alfven 
wave dispersion relations allows for a prediction of the 
observed resonance frequency. For the measured ex- 
perimental parameters, VA/wth,o ^ 1, suggesting that 
a kinetic calculation of the dispersion relation would 
be appropriate. However, because the collisionality 
is fairly high (Amfp,e ~ 0.2 m, fc||Amfp,c < 1) and 
in order to keep the calculation simple, a fluid dis- 
persion relation is used for kinetic Alfven waves |18j : 
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persion relation for the ion acoustic mode is Aw = 

Afc||Cs/yl -I- (Afcj^ps) . Three-wave matching relations 
predict that Aw = uj2 — wi, Afc|| = fc||2 -I- fc||i, and 
Afc_L = fc_L2 — fc_Li where 1 and 2 are subscripts asso- 
ciated with the counter-propagating Alfven waves. Some 
simple algebra and the assumption wi « a;2 = ^ 3> Aw 
lead to the equation: 
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FIG. 4. Resonance frequency A/ as a function of plasma and 
antenna parameters in both hydrogen and helium plasmas. 
Data is for a Langmuir probe a,t z = 4.47 cm in helium and 
z = 3.83 cm in hydrogen. The reference point is u = 480 kHz, 
Bo = 450 G in hydrogen and w = 220 kHz, Bo = 750 G in 
helium. k±ps is esimated from a Bessel function fit as 0.26 
in hydrogen and 0.28 in helium. Both magnetic field scans 
(circles) and antenna frequency scans (stars) are shown. The 
gray shaded region represents the value of A/ predicted by 
Eq. [l] taking into account uncertainties in the density and 
temperature measured by the Langmuir probe. 



Plugging in the experimental parameters used to produce 
Fig. |3J including a typical /3 ^ 4 x lO""*, Eq. [I] predicts a 
resonant frequency of 13 kHz. This agrees well with the 
experimental result. 

Eq. [l] is satisfied for a wide range of plasma and an- 
tenna parameters; this is shown in Fig. |4] For fixed ion 
mass, Eq.jT] implies that the resonant frequency is a func- 
tion of w/r^i. Therefore, magnetic field scans and scans of 
the main antenna frequency may be overplotted on Fig. [4] 
The resulting datapoints for each gas fall within the gray 
shaded region calculated using Eq. [ij The finite width of 
the gray region represents the statistical uncertainly from 
an average over similar Langmuir probe measurements. 
Temperature is measured by sweeping the voltage of a 
single tip; density is obtained from the measured tem- 
perature and ion saturation current measurements. 

To gain insight into the width of the response curve 
in Fig. I3J it is useful to model the non-linear interaction 
as a damped, driven oscillator. In the simplest possi- 
ble model, the ion acoustic perturbation is considered to 
be much smaller than the background density such that 
Pi I Pa ^ 1 and pi/po ^ vju). Self-consistent with this 
approximation, the parallel ion velocity perturbation is 
much smaller than the phase speed of the driven mode 
by the same order. Combining the MHD momentum and 
continuity equations and neglecting perpendicular prop- 
agation effects: 
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Eq. [2] describes a damped, driven oscillator system. 
The first and third terms represent the wave equation for 
an ion acoustic mode, the second term describes damping 
due to ion-neutral collisions, and the fourth term is the 
non-linear ponderomotive drive that results from interac- 
tion between the two Alfven waves. This acoustic mode 
drive term accelerates ions parallel to Bq through a non- 
linear vx B force in the parallel ion momentum equation. 
For the ordering assumptions used to derive Eq.[2]to hold, 
the amplitude of the Alfven wave drive at resonance must 
be small, b±/Bo <^ NyJjS, where N = v/ujq represents 
the collisionality normalized to the resonance frequency. 
During the time both Alfven waves are turned on, the 
system will respond at the drive frequency w/j = W2 — wi 



and the drive wave number fc^j 
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sponse function at w = ujd and k = ko follows from the 
linearization of Eq. [2J 
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where flo = i-^d/^o is the drive frequency normalized to 
the natural resonance of the system wq = koCg- Eq. Islis 
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FIG. 5. Spatial profile of the beat and Alfven waves in hy- 
drogen plasma. On the left plot, white vectors represent the 
superimposed magnetic fields from the two Alfven waves while 
the color scale shows spatial variations of the beat amplitude. 
The cathode and end mesh side antenneis are set to 480 kHz 
and 450 kHz respectively. Bo = 450 G. Magnetic data is av- 
eraged over the times at which the antenna-launched waves 
are in phase such that the magnetic field amplitude is at a 
maximum. The two current channels of the Alfven waves are 
visible in the upper right and lower left. In the right plot, 
the measured magnetic field amplitude and beat response are 
shown as a function of R, the distance from the orgin along 
the diagonal cut shown. The beat amplitude is peaked near 
the origin where the local magnetic field peaks. 



overplottcd on Fig. [s] for A^ = 0.14 (thin dash-dot line) 
and N ~ 0.36 (thin solid line). While the amplitude 
of the resonant response is well predicted by A^ = 0.14, 
the scaled inset figure shows that the best fit for the 
width of the peak is obtained for A^ = 0.36. Further- 
more, A^ = 0.36 is more than double the value A^ = 0.14 
obtained from the ringdown time in Fig. [2] This dis- 
crepancy is still under investigation and may be due to 
effects not included in the simple model, including finite 
perpendicular wave number and axial variation of plasma 
parameters. Despite this, the A^ = 0.14 result implies 
that the pondermotive force is of sufficient amplitude to 
drive the observed resonant response. 

The spatial profile of the beat response also suggests 
a pondermotive drive mechanism. This is shown exper- 
imentally in Fig. [5] The measured wave magnetic field 
vectors are plotted as white arrows; overlapping current 
channels for the two Alfven waves are indicated by the 
circulation pattern of these arrows in the upper left and 
lower right portions of the figure. The perpendicular 
wave number may be estimated by fitting the transverse 
magnetic field pattern of a single current channel to a 
spherical Bessel function of the first kind [Mj; for the 
plasma and antenna parameters associated with Fig. [5] 
this gives k± ^ 0.5 /cm and k±ps ~ 0.26. As indicated 



FIG. 6. Dispersion relation of beat waves at the resonance 
frequency. Experimental runs are the same as in Fig. |4] The 
best fit dashed line represents a phase speed of 29.1±0.7 km/s 
in the hydrogen plot and 12.5 ± 0.3 km/s in the helium plot. 



by the color scale, the beat amplitude is greatest near 
the origin which is where the Alfven wave magnetic field 
peaks. Qualitatively, this result agrees well with Eq. [3] 
Also consistent with the MHD theory presented, a scan 
of the antenna power (Fig. [2J Panel (D)) reveals that the 
beat-driven amplitude grows proportionally to the prod- 
uct of the two Alfven wave amplitudes. 

The resonance in the beat wave response is identified 
as an ion acoustic mode based on the dispersion relation. 
For each of the experimental runs in Fig. [4J the parallel 
wave number for the response at the resonance frequency 
is determined by examining the phase delay between two 
Langmuir probes closely spaced in z (One such experi- 
mental setup is shown in Fig. IT]). The result, shown in 
Fig. |6J is a linear dispersion relation with phase speed 
comparable to the sound speed for LAPD parameters. 
Using a kinetic dispersion relation [TT] for ion acoustic 
modes and assuming an ion temperature of 1 eV (pre- 
viously measured in both helium and argon |20j), the 
phase speed in helium requires Tg = 4.4 ± 0.3 eV, well 
in line with the value of 4.3 ± 0.9 eV obtained by ana- 
lyzing Langmuir probe sweeps during the beat wave time 
period. Applying the same set of assumptions to the hy- 
drogen data requires Tg = 6.7 ± 0.4 eV, well above the 
value of 4.3 ±1.0 eV measured. Since ion temperature 
has not been measured in hydrogen, one possible expla- 
nation for this discrepancy is T^ > 1 eV. This could be a 
combination of higher background ion temperature and 
enhanced ion heating by the launched Alfven waves in 
hydrogen plasmas, due to the lighter ion mass. For exam- 
ple, if the ion temperature is 2 eV, the ion acoustic mode 
dispersion relation [H] requires only Tg — 5.2 ± 0.3 eV 
which is within errorbar of the measured value. 

In summary, the first laboratory observations of the 



Alfvcn-acoustic mode coupling at the heart of parametric 
decay are presented. Counter-propagating Alfven waves 
launched from either end of the LAPD produce a reso- 
nant response identified as an ion acoustic mode based 
on the dispersion relation, spatial profile, and other fea- 
tures consistent with a simple MHD theory. Several areas 
for further investigation remain. Ion acoustic waves have 
never been directly launched by an antenna in the labora- 
tory at densities comparable to those in the LAPD. Thus, 
a new technique to directly launch ion acoustic waves is 
being developed and will be utilized for a detailed study 
of the damping mechanism. The new technique will also 
be used to investigate Alfven and ion acoustic wave cou- 
pling; the launched acoustic mode could potentially seed 
the parametric decay process. Additional studies may 
also focus on parametric decay from a single large ampli- 
tude Alfven wave; this is not possible under the present 
set of parameters due to insufficient Alfven wave ampli- 
tude. 
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